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Sur les solutions singulières des équations différentielles 

simultanées. 

Par B. Gouesat. 



La théorie des intégrales singulières des équations différentielles a donné 
lieu depuis Lagrange à un grand nombre de travaux ; l'un des plus connus est 
sans contredit le Mémoire de Mr. Darboux publié dans le Bulletin des Sciences 
Mathématiques (t. IV, 1873), dont les résultats sont aujourd'hui classiques. Je 
me suis proposé d'étendre les théorèmes de Mr. Darboux aux équations différen- 
tielles simultanées et aux équations d'ordre supérieur. La plupart des résultats 
auxquels on parvient ainsi pourraient, il est vrai, se déduire des propriétés établies 
dans le Mémoire du même auteur Sur les solutions singulières des équations aux 
dérivées partielles;* mais, comme ces résultats sont indépendants de la théorie 
des équations aux dérivées partielles, il m'a paru intéressant de les établir direc- 
tement. 

1. Pour ne pas interrompre la suite des raisonnements, je vais d'abord 
généraliser un des théorèmes fondamentaux sur les équations différentielles du 
premier ordre, démontré par Briot et Bouquet dans leur célèbre Mémoire sur ce 
sujet/}* Considérons un système d'équations différentielles de la forme 



Cil/ 

x ct = f( x ' y> z ) =ax + ty + cz + . 

dz 
*^-=/i (*,?/, a) = <hx+ % + C& + • 



(1) 



oùf(x, y, z) et fi{x, y, z) sont des fonctions holomorphes des variables ce, y, z 
dans le domaine du point x = y = z = 0. Proposons-nous de rechercher si ces 
équations admettent un système d'intégrales holomorphes s'évanouissant avec x . 

* Mémoires des Savants étrangers, t. XXVII. 

t Journal de V Ecole Polytechnique, XXXVI ème Cahier. 
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Si un tel système existe, il suffira, pour avoir les développements en séries, de 
calculer les dérivées successives de y et de z pour x = ; ces dérivées s'obtien- 
dront comme il suit. Après n différentiations successives, les équations (1) 

donnent d n+1 y d n y 3/ d n y df d n z i 

X dx n + 1 + n dx"~ dy daf + dz dx n + ' 

d n + h d»z_ _d£ d?y dj\ d^z 

X dx n+1 + n dx n ~ 'dy dx n + dz dx n + - " " 

si on y fait x = y = z = , on obtient les relations 
(n _ h) (*£\ - c (**-} = F \(ÈL\ (*£\ (*IÏL\ 1 ] 

" h (20 o + {n - Cl) (0)„= ^ Kir).' (S-)o' • • • • (aPQ J ' . 



(2) 



.F et i^ désignant des fonctions entières des dérivées 

/dy\ (d?y\ f d n ~h \ 

\dxJo' \dx*J '""\dx n - 1 Jo' 

/ d n v\ /d n z \ 
Ces relations déteïmineront ( -j\ ) et ( -=— ) au moyen des dérivées précédentes 

\dx /q \dix /o 

pourvu que le déterminant 

(n — b)(n — Ci) — bjp 

soit différent de zéro. Par conséquent, si l'équation 

(et — b)((ù — Ci) — èjc = 

n'admet pour racine aucun nombre entier positif, on pourra calculer de proche 
en proche toutes les dérivées successives 

/dy\ /<Py\ /dy\ 

\dxj ' \da?J ' \dx n J ' 

/ dz \ / dh \ / d n z \ 

\~dxJo' KdrfJo' ' ' ' "WVo ' 

de sorte que, si les équations (1) admettent un système d'intégrales holomorphes 
s'évanouissant avec x, ce système est unique et les intégrales seront représentées 
par les développements en séries 

y- i Kdxjo^ i.2\d«?j ^ — > 

x / dz \ x 2 / d?z \ 
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les coefficients étant calculés comme on vient de l'indiquer. Tout revient à 
démontrer la convergence de ces développements pour des valeurs de x de 
module suffisamment petit. Ici cette démonstration présente une difficulté 

spéciale, provenant de ce que les valeurs de (-j^j et (-rr) ne s'obtiennent pas 

par les seules opérations d'addition et de multiplication. En effet, en résolvant 
les équations (2), on aura en général dans les seconds membres des termes pré- 
cédés du signe — . On évite cette difficulté en opérant de la manière suivante. 
Ajoutons les équations proposées (1) après les avoir multipliées respectivement 
par deux constantes indéterminées % et fi; il vient 

x d ^ fXZ) - = (Jta + (ta 1 )x + iXb + (ib 1 )y + (Xc + (icjz + (3) 

Posons M> + ^ Xe + ^ 



■a, 



À /a 

ou X (6 — a) + (ibi = , 

%c-\- (ifa — o) = ; 

pour qu'on puisse satisfaire à ces relations par des valeurs de % et fi qui ne soient 
pas toutes nulles, il faudra que a soit racine de l'équation 

(o — b)(a — Cj) — b x c = . (4) 

Soit % une racine de cette équation ; on pourra satisfaire aux relations 

%{b — coi) + fibx = , 

"kG + [l ((?! (ùj) = 

en prenant pour X et (i des constantes dont l'Une au moins ne sera pas nulle. 
Supposons par exemple % différent de zéro ; si on pose ty + (tz = Y, on pourra 
prendre Y et z pour inconnues à la place de y et de z et le système (1) sera 
remplacé par le système 



dY 

x -=— = (%a + ^a a ) a; + o x F + ax % + . 



(5) 



œ i = a i œ + ~T F + ( c i~ &i -y) s + ai»* + j 

Remarquons encore que le coefficient de z dans la seconde équation est précisé- 
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ment la seconde racine % de l'équation (4). Si o 2 est différent de a t , on pourra 
poser , 

Z=-fY+(co,- ai )z 

et prendre Y et Z pour inconnues nouvelles ; alors les équations en Y et Z pren- 
drons la forme 

dY i _i_ v . 

d^ _ , , - , 
x -=— = a[x + 6) 2 Z + ...., 

mais cette seconde transformation, qui n'est pas toujours possible, n'est pas 
essentielle. Pour ne pas multiplier les notations, je suppose qu'on ait commencé 
par ramener le système primitif à la forme (5) ; alors le théorème général qu il 
s'agit de démontrer pourra s'énoncer ainsi. 
Les équations 



x -jL = ax + by + dx* + etf + fz % + . . . . =/(», y, z), 

x ~ = ajx + % + cjz + d&? + ejf + fz* + . . . . =/i(œ, y, z) , I 



(6) 



oùf(x, y, z),/ 1 (œ, 2/, 2) sowf efes fonctions kolomorphes dans h domaine du point 
x = y=z=0, admet un système d'intégrales holomorphes s' évanouissant avec x 
pourvu qu'aucun des coefficients b , c x ne soit égal a un nombre entier positif. 
En différentiant successivement les équations (6), on obtient 

x &y ■ fy _ df _■_ df dy 8/ ds 
da; 3 cfa; 3a; 8y da; 3s dx 

cPz dz _ 3/x 3/ t % 3/t dg 

dx 3 cfo da; 83/ da; 82 dx ' 

* da; 8 "^ da; 3 — da; 3 "^ " 1 ~ 83/ (fc 3 " r 3s da; 3 ' 

œ tf + tf~tf + "*" dy da; 3 "•" dz d3? ' 
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ce qui donne, en faisant x = y = z = , 

(•î).<i-»>=«. 



C de» A ( 2 &) ~" W A + " " * + 2 V32/3 Jo( ^ A V d* A ' 

( d?) (2 ~ Cl) = (a?), +•••• + 2 (a^) ( -£)H) + 6i ( dl) ' 



On en déduit 



1 \ dai /o 



/«y \ _ « / a^N _ 

U/o~ 1 — 6' UJo~ 

\âaïj— 2 — 6 

f<Pz\_ \d«?)j r ""~ i ~ KdydJoK dx J \ dx J^ ° l \ cM A 
\da?J Q — 2 — Cl 



(7) 



Pour démontrer la convergence des développements ainsi obtenus, il suffit 
d'employer un artifice tout-à-fait analogue à celui par lequel Briot et Bouquet 
ont démontré le théorème correspondant dans le cas d'une seule équation. Sup- 
posons que les fonctions f(x, y, z) et f x [x , y , z) soient développâmes en séries 
convergentes pour tous les systèmes de valeurs de x, y, z de module inférieur 
à r, et soit M une limite supérieure du module de ces fonctions dans ce domaine. 
Posons 

<p(x,u,v) = Ax + Bu + / _ A y i _^y i _ JL N -*{* + , } 



= Ax + Bu + M 



X 2 -f- W* + t> 3 -f" XU + XV 4- UV 



+ 



4>! (x, u, v) = A x x + BjU + Gjv + M | 



a; 8 -f- m 3 + a 3 -f- xu + aw + mc 

n>2 



+ . 



toutes des dérivées partielles des fonctions ty et «^ , à partir des secondes, sont 
réelles et positives pour x = m=v=0 et sont supérieures aux modules des dérivées 
correspondantes des fonctions / et f x pour x = ?/ = z = . Nous prendrons pour 
44 
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A, A lt B, les modules de a, %, b x respectivement et pour B, O x des nombres 
positifs inférieurs à l'unité qui seront déterminés plus loin d'une façon plus 
précise. 

Le système d'équations 

u — <p (ce , u , v) = 
v — $ 



> (ce , u , v) = ) 
i x (cc, w, v) = OJ 



(8) 



définit un système de fonctions holomorphes s'évanouissant avec ce, car le déter- 
minant fonctionnel des premiers membres se réduit à (1 — B)(l — 0,\ pour 
ce = u = v = . Dans un certain domaine autour de l'origine, ces fonctions seront 
représentées par des développements en séries convergentes 



u 



x / du \ a? / â?u \ , 

~~ T veto y fT2 V&V 



x / dv \ , af / et" 1 » \ 

v ~" T V &T/ 1T2 v &?/ 



(9) 



dont on pourra calculer tous les coefficients au moyen des équations (8). En 
différentiant ces relations on a, en effet, 

du 3<|> 3$ du dty dv 

dx ~~ dx du dx ' dv dx ' 



dv 3<£>i 3$! 

cfcc 3cc 3m c&c 

<fîu_(P<p 

dx % ~~ dx % "+"■••••• 

d 2 v 



du d<p x dv 

1 3v dx ' 



cZcc a 3x a ' 



+ 



+ 



3<|> #w 3$ ePv 

3« cfa; 2 ~ t ~ dv dx* ' 

3qt>! c^m s d$ x d 2 v 

3m cfcc 2 ' dv dx* ' 



On en déduit, en faisant ce = u = v = , 

cfe \ J. /_d»\ _ ^ 1 + ' Bl (^ / ) 

V e&c / o 1 — Ci 



tZw 



/ du\ 

V dx /q 1 — B 

,d?u\ _ WJ + • * • • + \dudv) Vdx-)y~fa) 
\d,x* ) — 



B 



\ (io) 



\dx* )~ 



_ V 3aVo + • ■ • • + 2 (g^gWoV «fa )„( dx ) + ^ ( 3cs ;i ) 



1-0, 
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On voit que tous ces coefficients sont réels et positifs. Puisqu'aucun des coeffi- 
cients b , c x n'est un nombre entier, le module de « — b , où n est un nombre 
entier positif, reste supérieur à une certaine limite m et, en prenant le nombre 
B tel que 1 — B soit inférieur à m, on aura constamment 

1_2*<|»_&|,.ou_^>|_JL ï 

on choisira de même le nombre positif C x de façon que Ton ait, pour toute valeur 
entière et positive de n, 

1 — G 1 <\n — c 1 \. 

Cela posé, la comparaison des formules (7) et (10) nous montre que l'on aura 



dy 
dx 

<$y 

dx 2 



o < \ dx X' 



.<( 

y I . / d'u \ I d*z I ^ / d" 
? o \~d¥)o' \tâ I J )\di 



du \ 
dx Jo 

<Pu\ 



dz 

dx 



d*v\ 



et, d'une manière générale, 



àTy 

dx' 



y\^(d n u\ I d n z I . /d n v\ 
" o \dp)o \dx» I o \dx")o 



Les séries (9) étant convergentes dans un certain domaine autour de l'origine, il 
en sera de même a fortiori des séries 



x fdy\ 



+ 



/ëy\ 



+...., 



1 KdxJo^ 1.2\cteVo 



x / dz \ z 3 / cPz\ 



ce qui démontre le théorème. 

Si nous revenons à la forme primitive des équations (1), nous pouvons 
énoncer la proposition générale suivante : 

Les équations (1) admettent un système d'intégrales holomorphes s 1 évanouissant 
avec x pourvu que V équation 

(a — b)(a — e^ — b x c =z 

n'admette pour racine aucun nombre entier positif. Un cas particulier de ce théo- 
rème a été démontré par Mr. Picard (Comptes-rendus, 1878). 
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(2). Nous allons maintenant étudier le cas où l'équation précédente admet 
pour racine un nombre entier positif; plusieurs cas sont à distinguer suivant 
qu'elle admet pour racines ou deux nombres entiers positifs. 

Supposons d'abord qu'il n'y ait qu'un seul nombre entier positif racine de 

l'équation (a — b)(a — Cj) — Z^c = ; 

soit p cette racine et soit o 2 la seconde racine. On pourra, comme on l'a vu plus 
haut, ramener le système d'équations proposées à la forme simple 

x -j- = ax + py + dx % + ...., 

dz 
x -r— = a x x + «23+ o\x % -\- .... 

On diminue les coefficients p et « 2 d'une unité en posant 

après p — 1 transformations de cette espèce, on sera ramené à un système d'équa- 
tions de la forme suivante 



x -p = ax + y + dx 2 + . . . 

x —r = a i x + G i z + dï 3 ? + • • • 
ax 



(H) 



où c x = Ma — p — 1 n'est pas un nombre entier positif. 

En général ce système n'admet pas d'intégrales holomorphes s'evanouissant 

avec x ; en effet, soit 

y = A ( pc + A 1 x i + 

z = B x + B 1 x i + 

les développements en séries de ces intégrales. En substituant dans les équa- 
tions (11), on aurait 

x (A + 2ÂjX +....) = ax + (A^ + A x x z +....) + œ 3 (. . . .) 
et par suite A = A + a ; 

ce qui exige que le coefficient a soit nul. Ainsi, lorsque le coefficient a est différent 
de zéro, les équations (11) n'admettent pas d'intégrales holomorphes s'evanouissant 
avec x. 
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Supposons maintenant que a soit nul ; si on pose 

y = hx, z = {ix, 
les équations (11) deviennent 

x-^=x{a + p% + y(i+ }, 

x -^ — aj + (i (Ci — 1) + x {a x + fta, + yifi + . . . . } ; 

la valeur initiale ^ de (i pour œ = doit vérifier la relation a x + /u ( c i — 1) = > 
tandis que la valeur initiale à de a peut être choisie arbitrairement. Si on 
prend pour \ une constante quelconque et si on pose 

% = \ + a/, ,« = r— ^ + ,«', 

on aboutit au nouveau système 

x^=zx{a> + pX + Y [i' + .. ..}, 

x^ = ( Cl - 1)^ + x{a{ + &*' + y^ + ....} 

auquel on peut appliquer le théorème général démontré plus haut. Ainsi, 
lorsque le coefficient a est nul, les équations (11) admettent ime infinité d'intégrales 
holomorphes s' évanouissant avec x; ces intégrales dépendent d'une constante arbi- 
traire îl . 

Supposons en second lieu que l'équation 

(a — b)(a — c 1 ) — bïC = 

admette pour racines deux nombres entiers positifs p et p + q (q > 0) . Par une 
série de transformations de même nature que les précédentes, on ramènera le 
système proposé à la forme 



x -p = ax + y + dx 2 +...., 

dz 
x -r~ = a x x + b x y + qz + d$? + 



(12) 



le coefficient Z> x pouvant être pris égal à zéro lorsque q est supérieur à l'unité. 
On démontre, comme tout-à-l'heure que ce système ne peut admettre d'intégrales 
holomorphes s'évanouissant avec x que si a est nul. 
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Si on a en même temps q=l, le système (12) devient 

x ^t=y + dx *+ » 

dz 
*> -fa = «i<k + % + » + dts& + 

Il peut arriver que l'on ait aussi b x = ; dans ce cas, en raisonnant comme plus 
haut, on démontrera que l'on doit avoir aussi a x = , et, si cette condition est 
satisfaite, on aura une infinité d'intégrales holomorphes s'évanouissant avec x car, 
en posant y=-%x, z = iix, on sera ramené à un système d'équations 

où les seconds membres sont des fonctions holomorphes de x, de à et de (i et on 
pourra choisir arbitrairement les valeurs initiales yl , (i des intégrales de ce 
système. 

Si b x n'est pas nul, le changement de variables 

y = hx, z = [ix 
conduit au nouveau système 

x-^=x(a + (& + yii+ ), 

x -J^- = a x + M + x (a x + /? x a, + y^ + ) . 

La valeur initiale /L de 2, sera donnée par la relation 

a x + b^ = , 
mais la valeur initiale (i de p reste arbitraire. En posant ensuite 

% — y + X', (i = ^ + (i', 

on arrive au système 

x^ = x(a' + P>X + y'ii> + ....), 

x ^ = b$ + x (ai + P& + y,V +•••■) 
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auquel on peut appliquer le théorème général. Il existe donc une infinité d'inté- 
grales holomorphes dépendant d'une constante arbitraire ^. 

Il reste à examiner le cas où q est un nombre entier supérieur à l'unité ; 
on peut se borner, comme on l'a. vu, à considérer le système 



dz 
x -r— = a x x + qz + dix 2 + 



(12') 



Si on pose y = x (/t + X) , z=zx (j-~Z. — *~ ^ ) ' ^° désignant une constante indé- 
terminée, on sera conduit à un système 

x -=- = x<p (x , /l , [i) , a; -^- = a a x + (ç — 1) (i + cZ 2 * 2 -f- • • • • 

Pour ce nouveau système, les deux racines de l'équation en a sont et g- — 1 ; si 
q est plus grand que 2, en diminuant successivement ces racines par le procédé 
qui a été expliqué plus haut de q — 2 unités on sera ramené finalement à un 
système de la forme 

x C ^~ = a M x+(2 — q)tt n) + , 

ax 

x^- = a{ n) x + !J. (n) + 

et ce système admettra des intégrales holomorphes s' évanouissant avec x pourvu 
que af ) soit nul. Comme ce coefficient dépendra en général de \, on voit qu'on 
aura une équation pour déterminer les valeurs convenables de % . Ayant pris 
pour \ une racine de cette équation, on aura une infinité d'intégrales holo- 
morphes dépendant d'une constante arbitraire j4 n) . 

3. Le théorème général qui vient d'être démontré s'étend sans difficulté à 
des systèmes formés d'un nombre quelconque d'équations de la forme suivante : 

X ~dV = aiX + hïVl + ^ + — + ^» + && + 

x ~~dV = a%x + htVx + c *y* + — + %» + P&? + ...., 

x ~ = a n x + h n y x + c n y, + + l n y n +p n x 2 + . . . . 
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Un pareil système d'équations admet toujours un système unique d'intégrales JioJo- 
morphes s' évanouissant avec x , pourvu qu'aucune des racines de l'équation 



\ — G) 



a 



h 



L—(ù 



= 



ne soit égale a un nombre entier positif. Si 0)1,03, . . . . a n sont les n racines de 
cette équation, on commencera par ramener le système proposé à la forme 

dyi 
dx 
dy-2. 



x-£ r =a l x + aji/i + . . . . + p 1 x 2 +...., 

a%x ,+ % x + o 2 y 2 + . . . . + pzx* + 



x 



dx 



x 



dy n 



dx 



= a n x + b n y x + c n y 2 + .... + a n y n + p n x % + . . . 



et la démonstration s'achèvera comme dans le cas de deux équations. On discu- 
tera de la même fapon le cas où l'équation en a admet pour racines un ou plu- 
sieurs nombres entiers positifs. 

On peut ramener à la forme précédente certains systèmes d'équations simul- 
tanées qui se présentent sous forme indéterminée. Considérons, pour fixer les 
idées, le système de trois équations 



dx 



dy 



dz 



zAxyz* ~ %A'x a yw ~ 2A"x a 'y'w 

où les trois dénominateurs sont nuls pour x = y = z=0. Pour avoir les inté- 
grales de ce système qui sont nulles pour x= 0, cherchons d'abord le degré infini- 
tésimal de y et de z par rapport à x . Soit % le degré de y et (i le degré de z ; 
supposons que les termes écrits au dénominateur soient les termes de moindre 
degré en x après qu'on a fait la substitution. Il faudra que l'on ait 

a + y,p + (iy = a' + %p' + uy' — % + 1 = a" + À/3" + py" — [i + 1 , 



ou 



a + *{(3 + 1) + [i (y + l) = a' + 1 +W + {i(y> + l) = a" + 1 + 1((3" + 1) + y". 
Imaginons que l'on ait marqué les points de l'espace de coordonnées 

(a, p+1, y+1), (o' + l, (3', y' + l), (a" + l, 0" + 1 , y"); 
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les relations précédentes expriment qu'il existe un plan 

X+XY+fiZ=p 

passant par trois de ces points au moins et laissant tous les autres points et 
l'origine de côtés différents. La détermination de ces plans s'effectuera sans 
difficulté par une construction géométrique tout-à-fait analogue à la construction 
plane bien connue pour développer les racines d'une équation algébrique. Sup- 
posons que l'on ait trouvé des valeurs de X et de fi répondant à la question ; ces 
valeurs seront commensurables et on pourra poser 

^ — JL „-JL 
r ' ™ r ' 

p, q, r étant trois nombres entiers positifs sans facteur commun. Si dans les 
équations proposées on fait le changement de variables 

cc = f, y = t p u, z = t g v, 

elles deviennent, après supression d'un facteur commun, 

du _ r (JW + ....) —pu (Au^v y + ....) 

dt ~ AuFv y + ' 

dv _ r (iW + ) — qv (Au^V +••••) 

* dt — Au^v y + . . . . 

Les valeurs initiales u , v de u et v pour t= sont déterminées par les équations 

simultanées 

r{A'u%vt + )— pu Q {Aulvl + ), 

r(A»uï"vl"+ ....)- qv,{Aulvl + ....). 

Soient u , v une solution de ce système ; en posant 

u = u + U, v = v + V, 

on sera ramené à un système de la forme (1), pourvu que l'on n'ait pas en même 
temps -4«o^o +....= 0; 

dans ce cas, on serait conduit à effectuer une nouvelle transformation de même 
espèce. On opérera de même dans le cas d'un nombre quelconque d'équations. 

4. L'étude des intégrales non holomorphes du système (1) serait sans doute 
très-intéressante ; mais, comme cette étude n'est pas essentielle pour la suite, je 
la laisse de côté dans ce travail. Je m'occuperai seulement du cas particulier où 
45 
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les racines de l'équation en a ont toutes les deux leur partie réelle négative. Il 
importe de définir exactement ce qu'on doit entendre par intégrales s'évanouis- 
sant avec ce. Je suppose que la variable ce tend vers l'origine suivant un chemin 
de longueur finie ayant à l'origine une tangente déterminée, de façon que l'argu- 
ment de x reste fini. En outre, je ne considère que des intégrales d'un degré 
infinitésimal déterminé par rapport à ce, c'est-à-dire qui peuvent se mettre sous 
la forme ce* (x + s) , % ayant sa partie réelle positive, x étant une constante diffé- 
rente de zéro et e une fonction de x qui tend vers zéro dans les mêmes conditions 
que la première. Avec ces restrictions, le système (1) n'admet pas d'autres intégrales 
s 1 évanouissant avec x que les intégrales Jiolomorphes lorsque les parties réelles des 
racines de l'équation en a sont négatives. 

Soient Ox, o 3 les deux racines de l'équation en o que je suppose distinctes 
pour fixer les idées ; le système (1) pourra se ramener à la forme 

x -~ = ax + a x y + dx* + ...., 
x -=— = a x x + (ù % z + dix % + . . . . 

Soient y — y x , z = «i les intégrales holomorphes de ce système. En posant 
y=zy x + u, z — z x + v et en remplaçant y 1 et z x par leurs développements en 
séries on est conduit aux équations 

x -7~--=.a l u + urf> (ce, u, v) + v4> {x, u, v), 

Q>V 

x — = a % v + ufa (ce , u , v) + v4>i (x, u, v) 

oùty, <£i, 4s ^î sont des séries ordonnées suivant les puissances croissantes de 
x, u, v sans termes constants. Je dis qu'un pareil système ne peut admettre 
d'intégrales de la forme 

u = x k (x + e) , v = af (x 1 + e') , 

les parties réelles de îl et de (i étant positives. En effet, supposons que la partie 
réelle de (j, soit égale ou supérieure à la partie réelle de %. En substituant les 
valeurs précédentes de u et de v dans la première des équations, on aura l'égalité 

de 

a,cc A (x + e)+cc A+1 -7- z=.iù 1 x x (x + e)-\-x x {x + e)q>(x, u, v) + x" (x 1 + e') i> (ce , u, v), 
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ou, en divisant par œ\ 

de 
(x + e)(à — %) = — x -z — \- (x + e)<p(x, u, v) + af~ A (V + e')4 , (», '«, «)• 

Imaginons maintenant que x tende vers zéro suivant un chemin de longueur finie 
ayant une tangente à l'origine ; les deux termes 

(x + e)<p(x, u, fc), x"~ K {x l + s')4'(x, w, v) 

de 
tendent vers zéro. Il doit en être de même du produit x -=— , au moins pour une 

certaine loi de décroissement de x. Supposons en effet que le module de 

de 
x -v- = n (x) reste constamment supérieur à une certaine limite m . De la relation 

de . . 

CC —r- = 7t (a?) 

dx v ' 

on tire dfÇ_ de 

x n (x) 

et if^L)^ r* lds 

\x / v x n (x) ' 

lorsque x tend vers zéro, le module du premier membre augmente indéfiniment 
tandis que le second conserverait une valeur finie ; ce qui est impossible. Le 
premier membre de la relation précédente doit donc tendre vers zéro avec x, 
ce qui ne peut avoir lieu que si % = «! . Mais alors la fonction ce* {x + s) ne 
tendrait plus vers zéro avec x puisque la partie réelle de à serait négative. On 
raisonnerait de même dans le cas où l'équation en a aurait ses racines égales. 

Il est à remarquer que le théorème n'est plus vrai si une seule des racines 
de l'équation en a a sa partie réelle négative. Par exemple, le système 
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admet l'intégrale 

II. 

5. Considérons une congruence de courbes, définie par les équations 

/(as, y, z, a, &) = 0,| ^ 

<p(x, y, z, a, 6) = J 
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où a et & désignent deux constantes arbitraires. Par chaque point de l'espace il 
passe en général un nombre fini de ces courbes ; car, si l'on remplace dans les 
équations (13) x, y, z par les coordonnées x , y , z de ce point, on a deux équa- 
tions pour déterminer a et b. Si l'on joint aux équations (13) la nouvelle équation 

df dcp df dcj> _ D{f,<p) _ (u) 

da ~db db da ~D(a, b)~~ ' K J 

on détermine sur chaque courbe de la congruence un certain nombre de points 
appelés points focaux ; le lieu de ces points focaux, quand on considère toutes les 
courbes de la congruence, forme en général une surface appelée surface focale, 
dont on obtiendra l'équation en éliminant a et b entre les équations (13) et (14) 
et qui se compose d'autant de nappes qu'il y a de points focaux sur chaque 
courbe de la congruence. Les courbes de la congruence sont tangentes en 
chacun de leurs points focaux aux différentes nappes de la surface focale. De 
plus, si l'on veut assembler les courbes de la congruence de façon à ce qu'elles 
aient une enveloppe, cette enveloppe sera située sur la surface focale et il est aisé 
de voir qu'elle sera définie par une équation différentielle du premier ordre. En 
effet, en chaque point A de la surface focale on connaît la direction de la tan- 
gente à la courbe de la congruence tangente à la surface en ce point et par suite 
la direction de la tangente à la courbe enveloppe qui passe par ce point.* 

Ke gardons maintenant, dans les équations (13), y et z comme des fonctions 
de la variable indépendante x définies par ces équations, a et b ayant des valeurs 
constantes quelconques. Une différentiation par rapport à x nous donne les deux 
nouvelles relations 

dx + dyV + dz z ~ ' 

dx dy y ' dz ' 

où , dy , dz 

a " dx ' dx 

En éliminant a et b entre ces équations et les équations (13) on est conduit à un 
système d'équations différentielles de la forme 



F(x,y,z,y',z') = 0, 
®(x, y, s, y', s') ; 



M (15) 

0;J 



* Pour la démonstration de ces propriétés, voir le tome II de la Théorie générale des surfaces de Mr. 
Darboux. Livre IV, Chapitre le'. 
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ces équations sont vérifiées par les fqnctions y et z de x définies par les équations 
(13), quelles que soient les valeurs constantes attribuées aux paramètres a et b. 
Nous dirons pour abréger que les courbes de la congruence sont les courbes 
intégrales du système (15), et forment l'intégrale générale. Mais ces équations 
(15) admettent en outre une infinité d'autres intégrales non comprises dans 
l'intégrale générale ; en effet, puisque ces équations établissent une relation entre 
un point de l'espace et la tangente à la courbe intégrale qui passe par ce point, 
il est clair que toute courbe tangente en chacun de ces points à une courbe inté- 
grale sera elle-même une courbe intégrale. Par suite, toutes les courbes enve- 
loppes des courbes de la congruence sont aussi des intégrales du système (15), et 
il est clair qu'en général elles ne font pas partie des courbes de la congruence ; 
nous les appellerons intégrales singulières. Nous sommes donc conduits à la 
proposition suivante : 

Les équations différentielles (15) admettent une infinité d'intégrales singulières, 
qui sont définies par une équation différentielle du premier ordre. 

6. Voici comment on pourra obtenir la surface focale et les solutions singu- 
lières en partant des équations différentielles elles-mêmes. Une des propriétés 
caractéristiques de la surface focale est la suivante ; par chaque point de cette 
surface passent deux courbes intégrales tangentes l'une à l'autre, une courbe de 
la congruence et une courbe enveloppe située sur la surface focale. Soient 
3%) 2/o> z o l es coordonnées d'un point de cette surface et y' , Zq les valeurs corre- 
spondantes de y' et de z/, relatives à ces deux intégrales. Le déterminant fonc- 
tionnel D ( J^ <£) 

devra être nul pour x=zx , y=y , z — z , y' = yd, z' — z[. En effet, si ce déter- 
minant était différent de zéro, on tirerait des équations (15) 

y'=:yl ) + P(x — x , y — y», z — z ), 
z' — zd + Q(x — x , y — y , z — zo), 

P et Q désignant des séries convergentes ordonnées suivant les puissances crois- 
santes de x — x , y — y , z — z , s'évanouissant pour œ = a3 , y = y , z = z . 
D'après le théorème fondamental de la théorie des équations différentielles, on 
en déduirait pour y et z des séries convergentes ordonnées suivant les puissances 
croissantes de x — cc , 

y = y + y' (x — x ) + a (x — x f + 

z = z + zl (x — a3 ) + a x (x — x f + . . . . ; 
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il n'y aurait donc qu'une seule courbe intégrale tangente au point (a: , y , z Q ) à 
la droite Z-Xp _ Y-y Q _ Z-z» . 

Par conséquent, il faudra que pour tout point de la surface focale les équations 
(15) et l'équation D (j. g) / 16) 

admettent une solution commune en y 1 , z 1 . Ainsi, on obtiendra l'équation de la 
sur/ace focale en éliminant y' et z 1 entre les équations (15)e£(16). 

Soit P (x, y, z) = l'équation ainsi obtenue ; les solutions singulières étant 
situées sur la surface focale, en éliminant z et 2/ entre les équations (15) et l'équa- 
tion de cette surface, on aura une certaine équation différentielle du premier 

ordre 4>{x,y,y') = 

qui définit les projections des solutions singulières sur le plan des xy. 

Si on étudie de plus près le procédé précédent, il donne lieu à plusieurs 
remarques. Des équations (15) on déduit, en différentiant, 



(17) 



dF , dF , , dF , , dF „ , dF .. n ■) 
9<ï> d® , . 3* . , 9* „ . d® _ 

-dx~ + -dy-y + te* + w f+ ^* Q] 

si on considère en particulier une intégrale singulière pour laquelle 

dFd^_dF d®_ 

dy' dzf dz 1 dy'~ 0, 

on pourra éliminer y" et z" entre les deux équations (17) et on obtient ainsi les 
deux nouvelles relations 

D(x, y') + D(y, y') & + D(z, y') z ~°> I 

D{F,* ) D{F,4>) D(F,*) \ K > 

D(x, z') + D{y, z 1 ) y """ D(a, a') a ~ U " J 

Les équations (16) et (18) se réduisent à deux équations distinctes; néanmoins, 
pour plus de symétrie, je les conserverai toutes les trois. On voit donc que pour 
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tout point de la surface P(x, y,z) — les équations (15), (16) et (18) doivent 
admettre une solution commune en y', z 1 . Or il est clair que, si les fonctions 
.F et <ï> sont quelconques, il n'existera pas de surface jouissant de cette propriété, 
et par conséquent il ne pourra exister d'intégrales singulières. 

On peut encore s'en rendre compte 'comme il suit. La première des équa- 
tions (15) exprime que la tangente à toute courbe intégrale passant en un point 
de coordonnées x, y, z est située sur un cône Payant pour équation 



F 



/ Y— y Z—z\_ n 



de même, la seconde équation (15) exprime que cette tangente est située sur un 
second cône T' ayant pour équation 

F- y Z 



^/ Y— y Z — z\_, 



En général, ces deux cônes T, T se coupent suivant un certain nombre de géné- 
ratrices distinctes et à chacune de ces génératrices correspond une courbe intégrale 
passant par le point de coordonnées x, y, z. Mais, si ce point vient sur la sur- 
face P(x,y,z) = Q, les deux cônes T et T' deviennent tangents suivant une 
génératrice commune G et, pour qu'il y eût une solution singulière passant par 
ce point, il faudrait évidemment que cette génératrice G fût située dans le plan 
tangent à la surface au point considéré, c'est-à-dire que l'on eût, pour tout point 

de cette surface, 

dP , dP . , dP , n , . 

te + -dy-y+!te z=:0 > < 19 ) 

y' et s' désignant les solutions communes des équations (15) et (16). Voici un 
moyen très-simple de s'assurer que cela n'a pas lieu en général. Remplaçons 
dans les équations (15) y' et z 1 par y' — m, z' — n, rn et n désignant deux con- 
stantes quelconques; nous ne changerons pas la surface P(x, y, z) et la relation 

(19) deviendra dP , dP , , s dP , , 

te + -dy-W- m ) + -dz~V- n ) = > 

relation qui ne pourra être vérifiée pour des valeurs quelconques de m et de w. 

Ainsi le procédé qui devrait conduire aux solutions singulières, appliqué a un 
système quelconque de la forme (15), ne fournit aucune solution. Un système quel- 
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conque d'équations différentielles simultanées n'admet pas d'une manière normale 
d'intégrales singulières. 

Ce théorème est tout-à-fait analogue, comme on voit, au théorème connu sur 
les équations différentielles du premier ordre,* et le paradoxe auquel on est con- 
duit s'explique de la même fapon. En effet, pour établir les propriétés des 
surfaces focales d'une congruence de courbes, on suppose implicitement, il est 
aisé de s'en assurer, que les fonctions / et <£ sont des fonctions continues de 
x, y, z, a, b dans le voisinage des valeurs des variables qui vérifient les relations 

Or, étant donné un système quelconque d'équations différentielles de la forme 
(15), on sait bien qu'il admet une infinité d'intégrales dépendant de deux con- 
stantes arbitraires ; mais rien ne prouve qu'on puisse mettre les équations des 
courbes intégrales sous la forme (13), les fonctions / et $ étant continues dans 
une étendue assez grande pour qu'on ait le droit d'appliquer la théorie des enve- 
loppes. Nous pouvons même affirmer, d'après ce qui précède, qu'il n'en sera pas 
ainsi en général. Ainsi, tandis que les systèmes d'équations simultanées formées 
directement par l'élimination des constantes admettent d'une manière normale 
une infinité d'intégrales singulières, au contraire un système d'équations différen- 
tielles pris arbitrairement ne pourra en avoir qu'à titre exceptionnel. Ceci nous 
montre qu'il est nécessaire, pour faire une théorie générale, de partir des équa- 
tions différentielles elles-mêmes et non de leurs intégrales. 

Il est bien clair d'ailleurs que les remarques précédentes n'enlèvent rien à 
l'intérêt des beaux théorèmes de Mr. Darboux sur les congruences de courbes, 
pas plus que nous ne devons rejeter la- théorie des enveloppes parce que les 
solutions d'une équation différentielle du premier ordre n'admettent pas en 
général de courbe enveloppe. 

7. Pour donner plus de précision aux. raisonnements, je supposerai que les 
fonctions F et <J> sont des fonctions algébriques entières et irréductibles de 
x, y, z, y', z' et, comme il ne saurait être question de passer en revue tous les cas 
particuliers qui peuvent se présenter, je n'examinerai que les hypothèses les plus 
générales. Je négligerai d'ailleurs toutes les difficultés provenant de valeurs 

* Bulletin des Sciences Mathématiques., t. IV, l êre Série, 1873. 
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infinies de y 1 et de z 1 , difficultés qu'on peut toujours faire disparaître par un 
changement de coordonnées. 

Pour un point quelconque de l'espace de coordonnées x , y , z , les deux 
cônes T, T qui ont pour équations 

Y— 2/ Z— v 






se coupent suivant p génératrices distinctes. Soient 

y— y» _ Z— 2b _ „ 

r - = — r— = JL — œ 

les équations d'une de ces génératrices ; le déterminant 

D{F, 3>) 

D(<y>,z>) 

n'étant pas nul pour x =x , y=y , z=2; , y , =yl ) , z , =z{ ) , les équations (15) peuvent 
être résolues par rapport à y 1 , z' et on en tire pour y 1 , z' des développements en 
séries convergentes ordonnées suivant les puissances de x — x , y — y , z — z , 

y' = yl + « (»— »o) + £ (2/ — 2/0) + r ( 2 — %) + — , 

rf = «J + ai (œ — flJ ) + ^i(«/ — 2/o) +n( 2 — z o) + 

On en déduit pour y et z les développements en série 

V = 2/o + y'o 0» — a- ) + h (a + £yj + ys$(aJ — £c ) a + . . . . , 
z = s + zj (x — x ) + h (%+ Piyô + yA)(% — x f+ . . . . 

qui seront convergents dans un certain domaine du point x . Ainsi, par un 
point quelconque de l'espace il passe en général p courbes intégrales, h tangentes 
distinctes, n'ayant en ce point aucune singularité. 

Il n'en est plus de même pour un point de coordonnées x , y , z pris sur la 
surface P(x, y, z) = dont on obtient l'équation en éliminant y' et 2/ entre les 
équations (15) et (16). En effet, pour tout point de cette surface, les deux cônes 
T, T deviennent tangents. Pour rester dans le cas le plus général, je supposerai 
que ces deux cônes sont simplement tangents le long d'une génératrice G et se 
coupent en outre suivant p — 2 génératrices distinctes. A. chacune de ces p — 2 
tangentes correspond- une courbe intégrale n'ayant au point considéré aucune 
46 
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singularité. Il nous reste à rechercher les courbes intégrales tangentes à la droite 
G au point cc , y , z . 

Supposons qu'on ait pris pour origine des coordonnées le point ce , y , z 
lui-même et la droite G pour axe des ce. Les équations (15) auront la forme 
suivante : A + Bf/ + CM + Dif + . . . 



oJ 



A l + B 1 y>+G l z' + D l yi°+... 

A, Ai, B, B lt . . . . désignant des polynômes entiers en ce, y, z tels que 
A = A 1 = BG 1 — GB 1 =0 pour x = y = z = 0. 

Les deux cônes T, T', relatifs à l'origine des coordonnées, auront pour équations 

(£)„—+ (0),-§- + .... = O, 

(B 1 ) ^+(G 1 \~+.... = 0, 

où (M) désigne, d'une manière générale, ce que devient M quand on y fait 
x = y = z = . Comme, par hypothèse, ces deux cônes sont simplement tan- 
gents, l'un au moins des coefficients (B) , (G) , (-Si)o, (Cï) sera différent de zéro. 
Supposons par exemple (B) #0; on pourra résoudre la première des équations 
(20) par rapport à y 1 et on en tirera pour y' une fonction holomorphe de ce, y, z, z' 
s'annulant en même temps que ces variables 

* = —%■* + .... 

et en portant cette valeur de y 1 dans la seconde des équations (20) on obtiendra 
une relation de la forme 

A + ( <?i - ^) * + M* + • • • • = o • 

Cette relation, considérée comme une équation en z', admet deux racines nulles et 
deux seulement pour ce = y = z = . Ces deux valeurs de z 1 peuvent être con- 
sidérées comme racines d'une équation du second degré dont les coefficients sont 
holomorphes en x, y, z dans le voisinage de l'origine. En résolvant cette équa- 
tion et en portant la valeur de z' dans l'expression de y', on obtiendra finalement 
pour y 1 et z 1 des expressions de la forme 



y' = A(x, y, z)+p(x, y, z)< /C(x, y, z) , 
d = B(x, y, z) + q(x, y, z)V G(x, y, z) 



;} (2i) 
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A(x,y,z), B(x,y,z), G(x,y,z), p(x,y,z), q(x,y,z) étant des fonctions 
holomorphes de x, y, z dans le voisinage de l'origine, et les trois premières 
s'évanouissant à l'origine ; nous admettrons de plus, ce qu'on peut toujours faire, 
que le radical doit être pris avec la même détermination dans les deux formules. 
Ces deux systèmes de valeurs de y 1 , z 1 deviennent égaux pour tout point dont les 
coordonnées vérifient la relation G(x, y, z)= 0; ceci nous montre que G(x, y, z) 
doit être identique à P(x,y,z) et, comme le plan tangent à cette surface 
P (x, y , z) = ne doit pas contenir l'axe des x, on aura 

C(x, y z) = ax -\-by + cz + dx 2 + . . . . 

le coefficient a étant différent de zéro. Soit 

A (x, y, z) = a^x + b^y + c x z + d^pt? + 

B (x, y, z) = a % x + b % y + c % z + d % v? +....; 

faisons le changement de variables 

x = x!\ y=- ux'*, z = va/ 3 , 
les équations (21) deviennent 



x 1 -j-j- = — 2m + 2a/ a (a x + b x u + CjW + . . .) + 2a;'p 1 (a/, u , v)*/a + bu + cv + 



dx' 
dv 
dat 



x 1 -T-j = — 2v + 2a3 ,a (a 2 + b 2 u + g s v ■+- . . .) + 2œ'q' 1 (a/, m, w)V f a + 6w -j- et? + • • • > 



les seconds membres de ces deux équations sont des fonctions holomorphes pour 
x' = u = v= Q puisque le coefficient a est essentiellement différent de zéro. 
D'après le théorème général, démontré au §1, ces équations admettent un système 
unique d'intégrales s'évanouissant avec x', représentées par les séries convergentes 

u = ax 1 + @x" + 

V — OLiX 1 + i^lSc' 2 + • • • • 

En revenant aux variables primitives x, y, z, nous voyons que les équations de 
la courbe intégrale sont les suivantes : 

x = f, 

y = afi + (3? + , 

Z = a i f + (3 1 t i + ; 

cette courbe gauche possède un point de rebroussement à l'origine. 
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Ainsi, la sur/ace P(x, y , z) = dont on obtient V équation en éliminant y 1 et z 1 
entre les trois équations 

F=0 *-0 ^'^ -Q 
Ji-U, cp_u, D y t j) -0 

est en général le lieu des points de rebroussement des courbes intégrales. 

En chaque point de cette surface passe une courbe intégrale ayant un rebrousse- 
ment en ce point, et la tangente de rebroussement est la génératrice de contact des deux 
cènes T, T' relatifs a ce point. 

8. Une intégrale des équations proposées sera dite singulière si pour chaque 
élément de cette intégrale on a à la fois 

l-V, <p_0, D y )S Q -o, 

de sorte que les théorèmes généraux sur les équations différentielles ne s'appli- 
quent plus à cette intégrale. Toute intégrale de cette espèce, s'il en existe, sera 
nécessairement située sur la surface P(x, y, z) = 0, obtenue en éliminant y 1 et z' 
entre les trois équations précédentes. Il faudra de plus que pour tout point de 
cette surface, au moins pour une certaine portion de cette surface, la génératrice 
de contact des deux cônes T, T, relatifs à ce point, soit située dans le plan tan- 
gent à la surface en ce point. Cette condition nécessaire est aussi suffisante ; 
car, si elle est remplie, on aura pour définir les solutions singulières une équation 
différentielle du premier ordre. Les fonctions F et 4> étant données, on pourra 
donc toujours reconnaître par des calculs algébriques s'il existe des solutions 
singulières et former l'équation différentielle qui les caractérise. 

Cette recherche pourra être facilitée dans certains cas par l'application de la 
règle suivante. Nous avons vu que pour tout point d'une solution singulière les 
équations (15), (16) et (18) devaient admettre une solution commune en 
y', z 1 . Réciproquement, supposons que pour tout point d'une certaine surface 
Q(x, y, z) = 0, (où Q (as, y, z) sera nécessairement un diviseur de P(x, y, z)) 
les équations 



D(F <ï>) 

F(x, y,z,m,n) = 0, <ï> (x, y, z, m, n) = 0, D x , ' / = 0, 

DjF,*) D(F,*) m , D(F t tb) 

D(x, m) ~*~ D(y, m) m + D(z, m) n - V> 

D(x,n) + D(y,n) m+ D(z,n) n ~ V 



(22) 
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admettent une solution commune en m, n. Quand on se déplace sur cette sur- 
face, m et n sont des fonctions de deux variables indépendantes, de x et de y par 
exemple. En différenciant totalement les premières équations (22), on obtient 

dF _ , dF . , dF . , dF 7 , 3F . 

-g — «as + -5— <*V H — g - ote + -g — dm + -g— an = , 
dx ' oy ff oz dm on 

D (F, <ï>) 

ou, en tenant compte de la relation ^w / = , 

' r D(m, n) ' 

i>(x, ot) dx + D(y, m) d y + D{z, m) dz ~ ' 

Admettons que l'un au moins des six déterminants fonctionnels qui figurent dans 

D (F, <ï>) 
ces formules, par exemple jyj r , ne soit pas nul ; le plan tangent à la surface 

Q (x, y, z) = aura pour équation 

et la quatrième des relations (22) exprime précisément que la droite 

X— x __ Y— y _ Z—z 
1 m n 

est située dans ce plan tangent ; et cette droite est la génératrice de contact des 

deux cônes T et T. On peut donc énoncer la proposition suivante : 

Lorsque pour tout point d'une surface Q (x, y, z) = les équations (22) admet- 
tent un système de solutions communes en m, n, sans que les six déterminants fonc- 
tionnels 

D(F,<b) D(F,<Î>) D (F, <£) B(F,^) D(F,<P) D(F,&) 
D(x, m) ' D(y, m) ' D(z, m) ' D (x, n) ' D(y, n) ' D(z, n) 

soient tous nuls pour ces valeurs de m et de n, les équations différentielles proposées 
admettent des solutions singulières. 
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Les intégrales singulières, si elles existent, peuvent être considérées comme 
les enveloppes d'autres intégrales. En d'autres termes, par chaque point de la 
surface Q(x, y, z) = il passe, outre une intégrale singulière, une intégrale non 
singulière tangente à la première. Imaginons que nous ayons pris pour origine 
des coordonnées un point de la surface Q(x, y, z) = et pour plan des xy le plan 
tangent en ce point, et soit 

«s = $(as, y) 

l'équation de cette surface. Si on change z en z-\-q>(x,y), il est clair qu'une 
transformation de cette nature n'altère pas les relations de contact entre deux 
courbes ; on peut donc supposer, pour la commodité du raisonnement, que la 
surface lieu des solutions singulières se réduit au plan des xy. Les deux systèmes 
de valeurs de y 1 , s' qui deviennent égaux pour z = seront donnés par des 
formules de la forme suivante 

y , = A(x,y,z)+p(x,y,z)*/z,) .^ 

z' = B(x, y, z) + q{x, y, z)Vz,) 

A(x, y, z), B (x, y , z) , p (x, y , z) , q(x, y, z) désignant des séries convergentes 
ordonnées suivant les puissances croissantes de x,y,z. En outre, comme z 1 
doit être nul pour tout point du plan des xy, B(x, y, z) doit contenir z en fac- 
teur dans tous ses termes et on peut écrire 

B{x, y, z) = zB l (x, y, z). 

Si dans les équations (23) on fait z = 0, elles se réduisent à l'équation unique 

^=A(x,y,0) 

qui détermine les solutions singulières. Si on pose d'autre part z = u 2 , elles 

deviennent dy .. „, . . „. 

-^= A(x,y,u*) + up(x,y,u*), 

^-^ = ^S l (x,y,u 2 )+ q(x,y,u 2 ), 

et dans le voisinage de tout point x = x , y=y§, u=0, ces équations admet- 
tent un système d'intégrales holomorphes pour lequel u n'est pas identiquement 
nul. Par chaque point (x , y ) du plan des xy passe donc une nouvelle courbe intê- 
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grale tangente h la solution singulière passant par ce point. Ce résultat est bien 
conforme à la théorie des solutions singulières faite en partant des congruences 
de courbes et de leurs surfaces focales. 

9. La proposition générale, qui fait le principal objet de ce travail, peut 
être établie sans avoir recours au théorème sur les équations différentielles qui 
a été démontré au début. Considérons d'abord une équation différentielle du 
premier ordre F(x,y, y 1 ) = (24) 

où F est un polynôme entier irréductible en x, y, y 1 ; intégrer cette équation, 
cela revient à exprimer x, y, y' eu fonction d'une seule variable indépendante de 
façon à vérifier l'équation (24) et la nouvelle relation 

dy — y'dx = . 

Or le système d'équations différentielles 

dx dy dy' 

■W-—W- ÔF ÏF , (25) 

dy 1 y dy' dx dy V 

admet l'intégrale première F(x, y, ?/) = const., et si l'on choisit les valeurs ini- 
tiales x , y , yô de fapon à ce qu'elles vérifient l'équation F(x , y , y' a )=- on 
aura aussi pour les intégrales 

F(x,y,y') = 0, dy — y'dx=0. 

L'intégration de l'équation (24) revient donc à l'intégration du système (25), les 
valeurs initiales vérifiant la relation 

F(x , y , y!,) = 0. 

Cela posé, supposons que pour x=-x , y=-yo l'équation (24) admette une racine 
multiple d'ordre n,y' = yç ) . Portons l'origine au point a; , y et prenons pour 
nouvel axe des x la droite de coefficient angulaire yô. On aura alors 

F(x, y, y') = P{x, y) + y ! P 1 (x, y) + y"P z (x, y) + . . . . +y"'P n (x, y) + . . . . 

les n polynômes P(x, y), P 1 (x, y), . . . . P n _ 1 (x, y) étant nuls pour a; = «/ = 0, 
et P n (x, y) n'étant pas nul pour x = y = 0. Les équations (25) deviennent 

dx 
P, {x, y) + .... + (n - l)2/'"- 2 P n _ 2 ( œ , y)+ .... 

dy dy* 

-tf\P 1 {x,y) + ....\-__d£_W ' 

dx dy t* 
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Plaçons-nous dans le cas général où il n'existe pas de solution singulière ; alors 

-p. h y' -p- ne sera pas nul pour x = y = y' = ; les équations précédentes 

pourront s'écrire 

—- = &(», y) + 2/Ç 2 ( œ , 3/) + . . . + ï/^Qn-tix, y)+y ,n ~ 1 Q n (x,y) + ..., 

^ r = y>Q 1 (x, y) + fQ,(x, y) + . . . + ^^-lO*. 2/) + f % Q« (* , 2/) + • • • . 

les coefficients Qi> ft» • • • • $«-i étant nuls pour œ = 2/ = o. Ces équations 
admettent un système d'intégrales holomorphes s'évanouissant avec y 1 , et il est 
aisé de voir que le développement de x commencera par un terme de degré n et 
celui de y par un terme de degré n + 1 . La courbe intégrale sera donc repré- 
sentée dans le voisinage de l'origine par les formules 

x = ay" 1 +••••, 
y = an/ n + 1 +.... 

Si n= 2, ce qui est le cas général, la courbe présente un point de rebroussement 
de première espèce; si ra> 2, l'origine est un point singulier d'ordre plus élevé. 
Etant donné un système d'équations simultanées 

F(x, y, z, y',d) = 0, ®(x, y, z, y',z') = Q, 

on verra comme tout-à-l'heure que l'intégration de ce système revient à celui du 
système 



dx dy dz 

D(F,d>) ~ , D(F,<Ï?) ~ d D(F,Q>) 
D&,*) y Dy,*) D{y\z') 

d£ d£ 

fl (ff*) ,,, »1ElÈÏ , :J D{F^) - D{F^) " Jpg .j DjF,®) 
~DV, x) + y Dy, y) ^ D(z>, z) D{x, y>) ^ D(y, y>) + « D(z, y') 



,(26) 



les valeurs initiales x , y , z , y^, z' Q vérifiant les relations 

F(x , y , z , yl, zQ = 0, $>(x , y , %, y' , z' a ) = 0. 

Supposons qu'en un point de coordonnées x , y , z les deux cônes T, T soient 
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tangents ; si on choisit ce point pour origine et la génératrice de contact pour 
axe des x, on aura 

F(x,y,z,t/,J) = A + Btf+ Œ + Dy»+ iEj/tt + HT + . . . . , 
$(x, y, z, y 1 , z') = a + by 1 + cz 1 + dy" + totf* +/a' 2 + . . . . 

ou A, B, G, D, E, F, . . . . a,b, c, d, e,f, . . . . sont des polynômes entiers en 
x , y , z , A et a étant nuls pour x = y = z=^0. Désignons d'une manière générale 
par L ce que devient une fonction quelconque L de x, y, z quand on y fait 
x = y = z = ; les deux cônes T, T' auront pour équations 

B °i + Ga Y + D * (x) + 2E °^ + F ° (t) + • • • • ' 

èo^ + co-^ + < *°( v xJ + ° x 5 "™vr/ + — 

Pour que ces deux cônes soient tangents suivant l'axe des x il faut que l'on ait 
B c — b G = 0. Nous pouvons même supposer que l'on a pris le plan tangent 
commun pour plan des xy. On aura alors 

B = & = , 

et nous admettrons que G et c ne sont pas nuls ; autrement l'axe des x serait 
une génératrice double pour l'un des cônes T, T'. Ecrivons les premiers termes 
des déterminants fonctionnels qui figurent dans les formules (26); 

D j}fj^] =Be— Gb + 2y'(Dc—Gd + Be — Eb) + %l{Eo— Ge+fB — bF) 

" d ''"' + 4 {De — Ed) f + 4 (D/— Fd) y'z 1 + 4 (Ef— Fe)z" + , 

+ <{•£-<'£+<-'»#} + •■•■•• 

D(F,Q) ,dA J«, , \,dB db dA on da\ 

+ n b dx-- B ^x- + 2e -dx-- 2E -dx-} + '--- 

T)( W <£i\ 
Pour x = y = z = 0, ^A ,' J se réduit à 

ly 1 (D c — G A) + 2 2' (-EoCo — <?cA>) + 

47 
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Le coefficient de y' ne sera nul que si on a D c — C d = 0, et cette relation 
exprime, il est aisé de s'en assurer, que les deux cônes T et T ont un contact du 
second ordre suivant l'axe des x. De même pour x = y = z = y' = z' = 0, les 
deux derniers dénominateurs des formules (26) se réduisent respectivement à 







dA 
dx 



dx /o 



et à zéro. Prenons le cas général où il n'existe pas de solution singulière ; alors, 
comme nous l'avons vu, on aura 



( 



dA _ da \ 

C dx-- G ^*° 



dxj 



et les équations (26) pourront s'écrire 



dx 

dy 1 



) = P X +Ptf + P s z> + 



dy __ 



dy' 



' T = P iy > + P % f + Pjy'd + 



dz 



dy 



7 = P 1 z> +Pty'z' +P s z" + 



(27) 



dz' 

ùf=Qi +Qœ' +Qs* +•••• 

P t et Qi étant des fonctions holomorphes de x, y, z dans le domaine de l'origine. 
Dans le cas général où les deux cônes T et T' sont simplement tangents suivant 
l'axe des x, P 1 et Q 1 sont nuls à l'origine mais P 2 est différent de zéro. Les 
équations (27) admettent un système d'intégrales holomorphes s'évanouissant 
avec y 1 et on aura pour les premiers termes du développement 

x = ay p + . . . . , 
y = Pf + 

z = yy" + , 

z> = èy" +.... 

Ces formules mettent en évidence le théorème général démontré plus haut (§7), 
mais nous voyons de plus que le plan tangent commun aux deux cônes T, T' 
rencontre la courbe intégrale en quatre points confondus avec l'origine. 
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Si on suppose maintenant que les deux cônes T et T aient un contact du 
second ordre suivant l'axe des x, P it P 2 et $1 seront nuls pour x=y=z= 0, 
et on voit facilement que les développements des intégrales en séries auront la 
forme suivante ™ — „,/• _i_ 

*/ = /%/<+...., 

z = wT + — . 
t = h" + ■■■■ 

La courbe intégrale présente à l'origine un point singulier d'espèce supérieure et 
le plan tangent commun aux deux cônes T, T' rencontre cette courbe en cinq 
points confondus à l'origine. Etant donné un système quelconque d'équations 
différentielles simultanées, les points singuliers précédents se présentent d'une 
manière normale. En effet, pour tout point de la surface P(x, y, z) = les 
deux cônes T et T' ont en général un contact du premier ordre. Si on veut que 
ces deux cônes aient un contact du second ordre, il faudra joindre à l'équation 
P(x, y, z) = une nouvelle équation de condition qui déterminera avec la pre- 
mière une courbe gauche pour tous les points de laquelle les deux cônes T et T' 
auront un contact du second ordre. 

On examinerait de même le cas où les deux cônes ont un contact d'ordre 
quelconque, ou le cas où l'un des cônes admet une génératrice double appar- 

tenant à l'autre et en général tous les cas où c -= G -^- n'est pas nul pour 

OX qX 

x = y = z=0. Si on avait en même temps (c -= (7-^-^ = 0, tous les 

\ d% dx J 

dénominateurs des équations (26) seraient nuls pour a; = ?/ = 3 = y = s' = 0, et 
on ne pourrait plus appliquer la méthode précédente. 

Il est facile de voir comment ces résultats se rattachent à la théorie des 
équations linéaires aux dérivées partielles. A chaque point de l'espace de coor- 
données (x, y, z) les équations (15) font correspondre m directions issues de ce 
point ; les surfaces telles que le plan tangent en chacun de leurs points contienne 
une des m droites correspondant à ce point seront définies par une équation aux 
dérivées partielles du premier ordre et du degré m en p, q, qui se décomposera en 
réalité en m équations linéaires. Les surfaces intégrales s'obtiennent, comme on 
sait, en associant suivant une loi arbitraire les courbes intégrales des équations 
(25). Si ces équations n'admettent pas de solutions singulières, comme c'est le 



360 Goubsat : Sur les solutions singulières des 

cas général, l'équation aux dérivées partielles n'aura pas non plus de solution 
singulière. Mais, si les équations (15) admettent des solutions singulières, la 
surface engendrée par ces courbes donnera une solution singulière de l'équation 
aux dérivées partielles. 
10. Exemple I. 

y — xy' = , œV 2 = x 2 + y* — 1 . 

Les deux valeurs de z' deviennent égales pour tous les points du cylindre 

à + y % — 1 = , 

y 

et la direction définie par la racine double est z 1 = , y' = — , c'est-à-dire la per- 

OC 

pendiculaire abaissée du point (x, y) sur l'axe des z. Cette direction n'est pas 
située dans le plan tangent ou cylindre, qui est par conséquent un lieu de points 
de rebroussement des courbes intégrales. Il est facile de le vérifier car l'inté- 
grale générale est 



y = O x x, z = Va 2 -f y % — 1 — arc tg */œ 2 + y % — 1 + G % . 

Exemple IL (V. Serret : Journal de Liowville, l ère série, t. XVIII, p. 29) . 

y = xy l + y' 1 + z 1 , z = z!x + y 1 *!. 

L'équation _> .' ,/ = est ici 
H D (y 1 , z 1 ) 

(x + y')(x+2y') — z' = 0. 

En éliminant y' et d entre ces trois équations, on trouve l'équation d'une surface 
du sixième ordre 

O 3 + Ixy — 9zf + 2 (x 2 + 3y)[l 2x (z — xy) — 2 (x 2 — y 2 ) 3 ] = . (28) 

Il est aisé de voir que, pour tout point de cette surface, les équations (15), (16) 
et (18) sont compatibles; nous nous trouvons donc dans le cas où il existe des 
solutions singulières. L'intégrale générale des équations proposées est en effet 
donnée par les équations 

y = ax + a 2 + b,) , 29 n 

z = bx + ab , ) 
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où a et b sont deux constantes arbitraires. Les courbes intégrales forment, 
comme on voit, un système de rayons rectilignes du second ordre et de la troisième 
classe; ces droites sont des tangentes doubles de la surface (28), qui contiendra 
par conséquent les arêtes de rebroussement des développables de la congruence. 
Pour obtenir ces développables, cherchons à déterminer b en fonction de a de 
fapon que les droites représentées par les équations (29) aient une enveloppe. 
On est conduit ainsi à l'équation différentielle 

/ db\* , db , „ ,„ . 

dont l'intégrale générale est 

b = Ga + G\ 

G désignant une constante arbitraire. L'arête de rebroussement correspondante 
aura pour équations 






l'équation (30) admet en outre l'intégrale singulière, 



»=-* 



et l'arête de rebroussement correspondante est 

y g-* 3 , 

_ ! 3 

Z -27" X - 



(31) 



(32) 



Les équations proposées admettent, outre l'intégrale générale (29), une infinité 
d'intégrales singulières représentées par les équations (31) et une intégrale singu- 
lière isolée (32) qui est l'enveloppe des premières. 

D'une manière générale, considérons le système d'équations simultanées 

F(y — xt/,z — sot, y', «0 = 0, * (y — »2/> a — a*'i î/, «0 = 0, 
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qui peut être considéré comme une généralisation de l'équation de Clairaut. On 
vérifie immédiatement que les équations (15), (16) et (18) se réduisent à trois 
équations distinctes ; le système doit donc admettre des intégrales singulières. 

En effet, l'intégrale générale est formée par les droites de la congruence 
F {y — ax, z — bx, a, b) = 0, <&(y — ax, z — bx, a, b) = 0, 

où a et b sont des paramètres arbitraires ; les arêtes de rebroussement des déve- 
loppables de la congruence seront les intégrales singulières. 



III. 

11. Les résultats qui précédent peuvent être appliqués aux équations différ- 
entielles du second ordre. Soit 

F(x,y,y',y") = (33) 

une équation du second ordre, où F est un polynôme entier irréductible en 

x, y, y 1 , y". Si on pose y' = z, l'équation (33) peut être remplacée par le 

système des deux équations 

F{x,y, z, O = 0,j , M) 

y' — z=0.) 
L'équation P(x, y, z) = , obtenue en éliminant z' entre les deux équations 

F(x,y,z,z') = 0, |^ = 

représente en général, nous l'avons vu, le lieu des points de rebroussement des 
courbes intégrales du système (34). Si on y remplace z par y' on a une certaine 
équation différentielle du premier ordre 

P(x,y,y') = 0, (35) 

que l'on obtiendrait évidemment en éliminant y" entre les deux équations 

F(x,y,y',y") = 0, J^ = 0, 

et dont il est aisé, d'après ce qui précède, d'avoir la signification. En effet, soient 
x, y, a les coordonnées d'un point M de la surface P(x, y , z) = 0; par ce point 
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passe une courbe intégrale du système (34) ayant un rebroussement en ce point 
et la droite X— x _ Y— y _Z— % 

~T~ y' ~~7~~ 

pour tangente de rebroussement. Cette courbe se projette sur le plan des xy 
suivant une courbe intégrale de l'équation (33) passant par le point (x, y), ayant 
un rebroussement en ce point et la droite de coefficient angulaire y 1 = z pour 
tangente de rebroussement. Par suite, en chaque point d'une courbe intégrale 
de l'équation (35) passe une courbe intégrale de l'équation (33) ayant un rebrous- 
sement en ce point et la tangente à la première courbe pour tangente de rebrous- 
sement. 

Ce théorème peut s'établir directement comme il suit. Supposons, ce 
qu'on peut toujours faire grâce à un changement de coordonnées, que pour 
x = y = y'= 0, l'équation (33) ait une racine double y" —-a. Pour des valeurs 
de x, y, y' suffisamment voisines de zéro, les deux valeurs de y" qui deviennent 
égales seront représentées par un développement de la forme 



y" = a + fax + \y + e#' +....) + Va^x + % + c 2 y' + . . . . ; (36) 

en posant, comme plus haut, y' = z, cette équation peut être remplacée par le 
système 



z' = a + fax + % -f <hz + ) + s/agc + % + c % z + ...., j 

y' = z. j 

Changeons encore z en ax + Z; le système devient 



Z = a x x -\-b^y + c x Z + . . . . + Va^x + h % y + c % Z +....,] 
y' = ax + Z. j 

A ce système on peut appliquer le théorème du n° 1 ; car, si on pose 

x = x'", y = x'^u , Z = x'*v , 



il devient 



X ŒU /a I #2 

— -r-y = U '+ aX. + X V, 

2 ax' 
-y ^-j = — v + x" fa + b x u + ci'» + . . . . \ : -f x' | a 2 + . . . .\ ; 
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ces équations admettent un système d'intégrales holomorphes s'évanouissant 
aveca; ' 2o, , 

v = T x + ' 

u = ~y x 1 ' + . . . . 

Par suite l'équation proposée admet une courbe intégrale tangente à l'axe des x 
à l'origine des coordonnées et représentée par les équations 

x = x 1 ', y = ~ô~ d* + ax 1 " + . . . . 

On voit que l'origine est un point de rebroussement de seconde espèce ; ce qu'on 
pouvait prévoir a priori puisque la dérivée seconde doit conserver en ce point 
une valeur finie. En réunissant tous ces résultats, on peut énoncer la proposi- 
tion ci-dessous : 

Etant donnée une équation différentielle du second ordre F(x, y, y 1 , y") = , en 

éliminant y" entre cette équation et la relation ^-jr = 0, on obtient une certaine équation 

différentielle du premier ordre P(x, y, y') = 0, dont les intégrales possèdent, en 
général, la propriété suivante. Par chaque point M d'une de ces courbes intégrales 
G il passe une courbe intégrale de Véquation F = , ayant un rebroussement de 
seconde espèce en M et la tangente en ce point a la courbe G pour tangente de 
rebroussement. 

On démontrerait aisément la proposition suivante, qui complète en quelque 
sorte la première : 

En égalant h zéro le coefficient de la plus haute puissance de y" dans Véquation 
F(x, y, y', 'y") = 0, on obtient une certaine équation différentielle du premier ordre 
Q(x, y, y') — dont les courbes intégrales sont telles qu'en général par un point 
quelconque M de Tune d 'elles G il passe une courbe intégrale de Véquation F = , 
ayant en M un rebroussement de première espèce et la tangente en ce point h la courbe 
G pour tangente de rebroussement. 

Si l'un des polynômes P(x,y, y 1 ), Q(x, y, y 1 ) ne contient pas y 1 , on devra 
regarder, suivant les idées de Clebsch, l'intégrale correspondante comme se com- 
posant d'un point quelconque de la courbe P(x, y) = ou Q (x, y) = et de 
toutes les droites passant par ce point. 
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12. Nous avons supposé jusqu'ici que le système auxiliaire (34) n'admet pas 
de solutions singulières. S'il en est autrement, toutes ces solutions singulières 
seront situées sur la surface P(x, y, a) = et leurs projections sur le plan des 
xy auront pour équation différentielle 

P(x,y,i/) = 0. 

Par chaque point de la surface P(x, y, z) = passent deux courbes intégrales 
du système (34) tangentes l'une à l'autre, dont les projections sur le plan des xy 
sont des courbes intégrales de l'équation (33) et, comme les valeurs de y' et de s' 
sont les mêmes pour les deux courbes dans l'espace, leurs projections sur le plan 
des xy auront un contact du second ordre. On voit donc que l'équation (33) 
admettra dans ce cas une infinité d'intégrales définies par une équation du 
premier ordre P(x, y, «/) = qui ont un contact du second ordre en chacun 
de leurs points avec une autre courbe intégrale. 

Les solutions singulières de l'équation (33) peuvent aussi être définies directe- 
ment, sans passer par l'intermédiaire du système simultané (34). Nous dirons 
qu'une intégrale de cette équation est singulière si pour tout point de cette 
intégrale la valeur correspondante de y" est racine multiple de l'équation 
F(x, y, y 1 , y") = 0. Il résulte de cette définition que de pareilles intégrales, si 
elles existent, devront vérifier l'équation P(x,y, y 1 ) = obtenue en éliminant 

O ET 

y" entre les équations F = et -^j, = . On sera ramené à rechercher s'il 

existe des intégrales communes aux deux équations 

F(x,y,y',y")=0, P{x,y,y')=0, 

ce qu'on pourra toujours faire par des calculs algébriques. Si l'équation F= 
est quelconque, ces deux équations n'auront pas d'intégrales communes et il 
n'existera pas de solutions singulières. Nous venons de voir quelle sera la pro- 
priété géométrique des courbes intégrales de l'équation P = . Mais, si les 
intégrales de l'équation P = appartiennent à l'équation F= , on pourra 
démontrer directement, comme plus haut, que ces courbes ont en chacun de 
leurs points un contact du second ordre avec une autre courbe intégrale. 

La recherche des solutions singulières peut être facilitée par l'application de 
la régie suivante. Il est clair d'abord que, s'il existe des solutions singulières, 
elles doivent satisfaire aux trois équations 
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F=o ^ = ^ + ^v + — v" = 0. 
* ' df ' dx ^ dy y ^ dy' y 



366 Goursat : Sur les solutions singulières des 

Inversement, si pour tous les points d'une courbe G ces trois équations admettent 
une solution commune en y", cette courbe est une solution singulière. En effet 
soit m la racine commune aux trois équations 

777/ I \ r, ^F ^ ^F i 3-F / , dF 

F(x,y,yJ, m) = 0,-^ = 0, _ + -_ y /+ _„, = ; 

quand on se déplace sur la courbe C, y, y 1 , m sont des fonctions de x dont les 
dérivées vérifient la relation suivante 

i^ _|_ i^! '4-i^ // i dF dm _ 
dx i ~ dy y dy 1 y "*" dm dx ' 

qui, comparée avec les premières, se réduit à 

On aura donc y" = m, à moins que l'on n'ait aussi -~-j = 0. 

Remarque. Si l'équation .F= admet des solutions singulières définies par 
l'équation du premier ordre P = , cette équation P = pourra avoir elle-même 
une solution singulière et cette solution singulière n'appartiendra pas en général 
à l'équation ¥■=■ . Car cette solution n'a en général qu'un contact du premier 
ordre avec les autres intégrales de P = . 

13. Exemple I. y"' — 2y'y" +1 = 0. 

L'équation P (x, y, y') = est ici y'* — 1 = 0; l'intégrale générale se compose 
de lignes droites y=±x+C 

et il est facile de vérifier que ces lignes droites sont des lieux de points de 
rebroussement de seconde espèce pour les courbes intégrales de l'équation pro- 
posée ; l'intégrale générale, que l'on trouve aisément, est en effet représentée par 
les équations 



p~ — 1 

X = 2 " 


-— pVp* — 1 


+ 


f*< 


p 8 — 1 
V= 3 " 


(p 2 — 1) \Zp a - 
3 


-1 


+ 0,. 


Exemple II. 









(1 + x *) y »>- ( 2 xy' + ~)y" + y" + xy>-y= .* (37) 

* Lagrange, Leçons sur le calcul des fonctions. Serret, Cours d'Analyse, t. II, p. 895. 
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Les deux valeurs de y" deviennent égales si on a 

On vérifie sans difficulté que les intégrales de l'équation (38) vérifient l'équation 

(37) ; il y a donc des solutions singulières. L'intégrale générale de l'équation 

(37) est en effet 

y =r an? + bx + 4a 3 + b % 

et celle de l'équation (38) 



V 1 Qy + 4îc 2 + x i — x «/ 1 + a? — log (x + V 1 + as 2 ) = H. 
L'équation (38) admet en outre une intégrale singulière 

X % X* 

y— 4 -yT 

qui n'appartient pas à l'équation (37). 

14. Il est aisé d'étendre ces considérations à des systèmes d'un nombre quel- 
conque d'équations du premier ordre tels que 



-^i(œ» yi, y% yn\ y'i, yi yL) = o» 

-*»(»» y u y t , y n ; y{, y' % y' n ) = 0, 

F n {x, yn y* yn) yi yi, — yL) = o. . 



(39) 



Appelons point dans l'espace à (n + 1) dimensions tout système de valeurs par- 
ticulières x Q , yl, . . . . yl attribuées aux variables x, y it courbe l'ensemble des 
points dont les coordonnées sont fonctions d'une seule variable indépendante, 
sur/ace l'ensemble des points qui vérifient une seule relation 

*(»i yi> ya 2/») = 0. 

Tant que le déterminant fonctionnel 

D{F 1 ,F % F n ) 

D(y'\, yi, — yL) 

sera différent de zéro, on pourra résoudre les équations (39) par rapport à 
y{> yi> • • • • y'n e * appliquer le théorème fondamental de Oauchy. Mais il n'en 
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sera plus de même si ce déterminant est nul. En éliminant y{ y n entre les 

équations (39) et l'équation 



DJF^F, F n ) _ 



0, (40) 



D(y[,yi, y'n) 

on obtient une certaine surface 

P{x,yi,y*, Vn) = 0. (41) 

En général, par chaque point de cette surface passe une courbe intégrale pré- 
sentant -en ce point une singularité que l'on peut définir comme il suit ; par une 
substitution linéaire convenable, x, y lt . . . . y n sont représentées par des déve- 
loppements de la forme 

x — t\ 

y x = a/ -h 

y % — a$ + 



y n = a n t»+*+.... 

Si, en chaque point de cette surface, la direction particulière définie par les équa- 
tions (39) et (40) est située dans la variété linéaire à n dimensions tangente à cette 
surface, les conclusions sont tout-à-fait différentes. La surface P = est un 
lieu d'intégrales singulières, que l'on peut définir par un système de {n — 1) 
équations différentielles du premier ordre ou, ce qui revient au même, par une 
équation différentielle unique d'ordre (n — 1). 

Pour appliquer cette théorie à un exemple, reprenons un des problèmes 
traités par Serret ;* proposons-nous de déterminer une courbe gauche, connaissant 
la courbe lieu des centres de courbure. Si on choisit x comme variable indé- 
pendante, on aura à rechercher les- intégrales de deux équations différentielles du 
second ordre 

F(x,y, z,y', z>, y", z") = 0,\ (42) 

<ï>(a5, y, z, y', z', y", z") = 0;i 

* Journal de IÀouville, l ère Série, t. XVIII, p. 23. 
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ce système peut être remplacé par le suivant, qui ne contient que les dérivées du 
premier ordre 

F(x, y, z, u, v, u', «0 = 0, 
Q>(x, y, z, u, v, u', v') = 0) 
y' — «= 0, 
z 1 — v= 0. 



} (43) 



L'intégrale générale doit contenir quatre constantes arbitraires ; or on connaît 
immédiatement des courbes intégrales répondant à la question, les cercles qui 
ont leur centre en un point quelconque de la courbo donnée G. En dehors de 
ces solutions, évidentes a priori, il existe d'autres courbes répondant à la question 
que l'on pourra définir directement. Considérons une développable passant par 
la courbe G; cette développable sera définie si on se donne l'angle ty que fait en 
un point quelconque M de G le plan tangent à cette développable avec le plan 
osculateur à la courbe G en ce point. Soit D une développable obtenue ainsi, 
et G la génératrice passant en M. Toute courbe gauche répondant à la question 
peut évidemment être considérée comme l'arête de rebroussement de la surface 
enveloppe d'un plan P mené par M perpendiculairement à la génératrice G d'une 
certaine développable D; et il faudra de plus que le point de contact du plan 
P avec son enveloppe soit situé dans le plan tangent à la développable D suivant 
la génératrice G. En exprimant cette propriété, on est conduit, il est facile de 
s'en assurer, à une équation différentielle du troisième ordre dont les intégrales 
fournissent la véritable solution du problème proposé. Ces courbes correspondent 
à des solutions singulières du système (43) ; en effet, par chaque point de l'une 
d'elles passe une courbe intégrale ayant avec elle un contact du second ordre, le 
cercle osculateur lui-même, et pour ces deux intégrales x, y, z, u, v, y', z 1 , u 1 , v 1 
auront les mêmes valeurs. 

On peut encore s'en rendre compte autrement. Cherchons les courbes inté- 
grales passant par un point donné M de l'espace et tangentes à une droite donnée 
MM ' passant par ce point. Autrement dit, cherchons les intégrales du système 
(43) qui correspondent à des valeurs initiales données de x, y, z, u, v. Pour 
cela, au point M nous mènerons le plan P perpendiculaire à la droite MM 1 , et 
nous prendrons les points d'intersection Q 1 , 2 , .... 0p de ce plan avec la courbe 
G; si ces p points sont distincts, ce qui est le cas général, les seules courbes 
répondant à la question seront les cercles décrits des points X , 8 , . . . . 0^ comme 
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centres et tangents à la droite MM'. Pour qu'il y ait une véritable courbe inté- 
grale tangente en M à la droite MM', il faudra que le plan P soit tangent à la 
courbe (7; par suite deux des systèmes de valeurs de u', v', fournies par les équa- 
tions (43) seront venu se confondre : ce qui est la propriété caractéristique des 
solutions singulières. 

15. Prenons encore une équation unique d'ordre n 

F(x,y,y',y",....y™)=:0, (44) 

ou F désigne une fonction entière de x, y, y', . . . . y in) ; si, pour 

x = x , y = y y n ~ 1 z=y^- 1 , 

l'équation (44) admet une racine multiple y" = y%, on ne peut plus appliquer les 
théorèmes généraux. Pour trouver la singularité correspondante de l'intégrale, 
imaginons, ce qu'on peut toujours faire, que cette circonstance se présente pour 
les valeurs initiales cc = y Q = yô = ; et supposons d'abord que l'on ait aussi 
y' ' = yl" = ....= y* = o , la valeur de y {n) qui se réduit à zéro étant racine 
double seulement. Les deux valeurs de y (n) qui deviennent nulles seront repré- 
sentées par un développement de la forme 



y [n) = «i» + a % y + «32/ + . . . + a n+J y {n - 1} + . . . + Vb x x + \y +%' + ... (45) 
Posons y' = %, y" = u% , . . . . j/ (n-1) = u n _- L ; 

on peut remplacer l'équation (45) par le système 



«1 = «2, 



u n— 2 — u n — 1) 



<_i = ajx + <hy + a 3 «i + + ^hx + % + è 3 % + . . . 

Faisons encore le changement de variables 

x=xf\ y = x"Y, u- l = x'*Ui, . . . . u n _ 1 = vl^Un-i) 



(46) 



équations différentielles simultanées. 371 

le système (46) devient 



4f + r=«^. 



x' dU„_ 2 



2 &/ 



(47) 



D'après le théorème général du n° 1, ce système admet un système d'intégrales 
holomorphes s'évanouissant avec x', et on vérifie aisément que le développement 
de F commencera par un terme en (x') 2n + 1 . On aura donc pour y un développe- 
ment suivant les puissances de x i de la forme suivante 

y = ax n+i + (3x n + 1 + iyx n + *+ .... 

Supposons maintenant que pour x = y = y' z= o, y" = yô' . . . . y (n ~ 1] = yk n ~ l \ 
l'équation (44) admette une racine double y {n) = y { ^\ On ramènera ce cas au pré- 
cédent en posant 

y— 2 œ +i. 2 .3 a!+ +1.2 n x +z > 

on en conclut que le développement de y sera de la forme 

y = \ x* + ïf^g x 3 + . . . . + 1-2 .. ..„ «" + «a n + * + /fo" +1 + • • • • (48) 

La courbe intégrale présente à l'origine un point singulier d'espèce particulière, 
analogue à un point de rebroussement de seconde espèce. 

Pour que l'équation (44) ait une racine double en y (n \ il faut que l'on ait en 

même temps ^p 



dy 



■(») 



= 0; 



l'élimination de y in) entre ces deux équations conduit à une relation entre y et 
ses dérivées jusqu'à l'ordre n — 1 . 

*(«,?, y' y^~ 1 >) = 0. (49) 
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D'après ce que nous venons de voir, les courbes intégrales de cette équation (49) 
possèdent la propriété suivante. En chaque point M de l'une d'elles G, il passe 
une courbe intégrale de V équation proposée (44) ayant avec la courbe G un contact 
d) ordre n — 1 et présentant en ce point un point singulier de V espèce caractérisée par 
le développement (48). 

Si les intégrales de l'équation (49) vérifient aussi l'équation (44), ce qui ne 
peut avoir lieu que dans des cas particuliers, ces intégrales seront des intégrales 
singulières. On démontrerait comme précédemment qu'en chaque point de l'une 
d'elles passe une seconde courbe intégrale ayant un contact d'ordre n avec 
l'intégrale singulière. 

Pakis, Janvier 1889. 



